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Abstract 

A search for excited quarks is performed using the full e ± p data sample collected by the 
HI experiment at HERA, corresponding to a total integrated luminosity of 475 pb _1 . The 
electroweak decays of excited quarks q* — » qj, q* — > qZ and q*-^qW with subsequent 
hadronic or leptonic decays of the W and Z bosons are considered. No evidence for first 
generation excited quark production is found. Mass dependent exclusion limits on q* pro- 
duction cross sections and on the ratio //A of the coupling to the compositeness scale are 
derived within gauge mediated models. These limits extend the excluded region compared 
to previous excited quark searches. 
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1 Introduction 



The three-family structure and mass hierarchy of the known fermions is one of the most puzzling 
characteristics of the Standard Model (SM) of particle physics. Attractive explanations are 
provided by models assuming a composite structure of quarks and leptons [1]. The existence of 
excited states of leptons and quarks is a natural consequence of these models and their discovery 
would be a convincing evidence of a new scale of matter. Electron ^proton interactions at high 
energies provide the opportunity to search for excited states of first generation fermions. For 
instance, excited quarks (q*) could be singly produced through the exchange of a 7 or a Z boson. 

In this paper a search for excited quarks using the complete e ± p HERA collider data of the 
HI experiment is presented. Electroweak decays of the excited quark into a SM quark and a 
gauge boson (7, W and Z) are considered and both hadronic and leptonic decays of the W and 
Z bosons are analysed. 

The data were recorded at an electron beam energy of 27.6 GeV and proton beam energies 
of 820 GeV or 920 GeV, corresponding to centre-of-mass energies y/s of 301 GeV and 319 GeV, 
respectively. The total integrated luminosity of the data sample is 475 pb _1 . The data comprise 
184 pb _1 recorded in e~p collisions and 291 pb _1 in e + p collisions, of which 35 pb _1 were 
recorded at \fs = 301 GeV. With more than a twelve-fold increase in statistics and a higher 
centre-of-mass energy, this analysis supersedes the result of previous searches for excited quarks 
at HERA by the HI [2] and ZEUS [3] Collaborations and is complementary to exclusion limits 
obtained at the LEP collider [4] and at the Tevatron [5-8]. The analysis also complements 
searches for first generation excited neutrinos [9] and electrons [10] at HERA. 



2 Phenomenology 

In the present study a model [11-13] is considered in which excited fermions are assumed to 
have spin 1/2 and weak isospin 1/2. The left-handed and right-handed components of the 
excited fermions form weak iso-doublets F£ and F^. Interactions between excited and ordinary 
fermions may be mediated by gauge bosons, as described by the effective Lagrangian [12, 13]: 



'int. ~ 2A - t R 



F* n o» v 



F L + h.c. . (1) 



Only the right-handed component of the excited fermions F^ is allowed to couple to light 
fermions F L , in order to protect the light leptons from radiatively acquiring a large anomalous 
magnetic moment [14, 15]. The matrix a^ v is the covariant bilinear tensor, W^ v , B^ v and G^ u 
are the field- strength tensors of the SU(2), U(l) and SU(3)c gauge fields, r a , Y and A a are 
the Pauli matrices, the weak hypercharge operator and the Gell-Mann matrices. The standard 
electroweak and strong gauge couplings are denoted by g, g' and g s , respectively. The parameter 
A has units of energy and can be regarded as the compositeness scale which reflects the range 



'in this paper the term "electron" is used generically to refer to both electrons and positrons, unless otherwise 
stated. 
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of the new confinement force. The constants /, /' and f s are. coupling parameters associated 
with the three gauge groups and are determined by the yet unknown composite dynamics. 

Following this model of gauge mediated interactions, excited quarks may be resonantly pro- 
duced in ep collisions through a gauge boson exchange between the incoming electron and a 
quark of the proton (see figure 1(a)). The exchange of excited quarks in the w-channel (fig- 
ure 1(b)) is also possible for the high q* masses and low values of A probed in this analysis. 
For example, for A = 50 GeV, f — f — 1 and an excited quark of mass M q * = 300 GeV, the 
w-channel production cross section is 0.016 pb while the resonant production cross section is 
0.27 pb. 




Figure 1: Diagrams for the production and decay of excited quarks in ep collisions. 

The excited quark may decay into a quark and a gauge boson via q*— >g7, q*^qW, q*^qZ 
and q*^qg. For a given q* mass value and assuming a numerical relation between /, /' and f s , 
the q* branching ratios are fixed and the production cross section depends only on //A. Only 7, 
W and Z decays of the q* are considered in the present study. In this analysis, the assumptions 
are made that the coupling parameters / and /' are of comparable strength, with the relationship 
/ = /', and that f s = 0. These assumptions lead to results which are complementary to the q* 
searches performed at the Tevatron [5-8], since at a pp collider excited quarks are dominantly 
produced in a quark-gluon fusion mechanism, which requires f s 7^ 0. The effect of non-zero 
values of f s is also studied in the present analysis. 



3 Simulation of Signal and Background Processes 

A Monte Carlo (MC) program developed for this analysis is used for the calculation of the q* 
production cross section and the simulation of signal events. The events are simulated using 
the cross section calculated from the Lagrangian described in equation (1) using the CompHEP 
program [16]. Both resonant q* production and u-channel exchange processes, as well as their 
interference are included. Initial state radiation of a photon from the incident electron is in- 
cluded using the Weizsacker- Williams approximation [17]. The proton parton densities are 
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taken from the CTEQ5L [18] parametrisation and are evaluated at the scale = y/sx, where 
x is the momentum fraction of the proton carried by the interacting quark. The parton shower 
approach [19] is applied in order to simulate Quantum Chromodynamics (QCD) corrections 
in the initial and final states. Hadronisation is performed using Lund string fragmentation as 
implemented in PYTHIA [19]. 

In the MC generator the full transition matrix including both q* production and decay is 
implemented. This is important if the natural width of the q* is large, which is typically the case 
at high mass. In order to incorporate the width effects in the signal efficiency determination, 
events are generated with the coupling //A corresponding, for each q* mass, to the expected 
boundary of the probed domain in the plane defined by M q * and //A. 

Excited quarks will be searched for in the 57, qqq, qeu, q/iu, qee and q/i/i final states. The 
SM background processes that may mimic the q* signal are prompt photon production, neutral 
current (NC) deep inelastic scattering (DIS), photoproduction, single W boson production and 
lepton pair production. 

The RAPGAP [20] event generator, which implements the Born, QCD Compton and Boson 
Gluon Fusion matrix elements, is used to model NC DIS events. The QED radiative effects 
arising from real photon emission from both the incoming and outgoing electrons are simulated 
using the HERACLES [21] program. Direct and resolved photoproduction of jets and prompt 
photon production are simulated using the PYTHIA event generator. The simulation is based 
on Born level scattering matrix elements with radiative QED corrections. In RAPGAP and 
PYTHIA, jet production from higher order QCD radiation is simulated using leading logarith- 
mic parton showers and hadronisation is modelled with Lund string fragmentation. The leading 
order MC prediction of NC DIS and photoproduction processes with two or more high trans- 
verse momentum jets is scaled by a factor of 1.2 to account for the incomplete description of 
higher orders in the MC generators [22,23]. Charged current DIS events are simulated using the 
DJANGO [24] program, which includes first order leptonic QED radiative corrections based on 
HERACLES. The production of two or more jets in DJANGO is accounted for using the colour 
dipole model [25]. Contributions from elastic and quasi-elastic QED Compton scattering are 
simulated with the WABGEN [26] generator. Contributions arising from the production of sin- 
gle W bosons and multi-lepton events are modelled using the EPVEC [27] and GRAPE [28] 
event generators, respectively. 

Generated events are passed through a GEANT [29] based simulation of the HI apparatus, 
which takes into account the actual running conditions of the data taking, and are reconstructed 
and analysed using the same program chain as is used for the data. 



4 Experimental Conditions 

A detailed description of the HI detector can be found in [30]. Only the detector components 
relevant to the present analysis are briefly described here. The origin of the HI coordinate 
system is the nominal ep interaction point, with the direction of the proton beam defining the 
positive z-axis (forward region). Transverse momentum (P T ) is measured in the x-y plane. The 
pseudorapidity r\ is related to the polar angle 9 by rj = — In tan(#/2). The Liquid Argon (LAr) 
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calorimeter [31] is used to measure energy and direction of electrons, photons and hadrons. 
It covers the polar angle range 4° < 9 < 154° with full azimuthal acceptance. Electromag- 
netic shower energies are measured with a precision of a(E)/E = 12%/v/E/GeV © 1% and 
hadronic energies witha(E)/E = 50%/v/E/GeV© 2%, as measured in test beams [32,33]. In 
the backward region, energy measurements are provided by a lead/scintillating-fibre (SpaCal) 
calorimeter 2 [34] covering the angular range 155° < 9 < 178°. The central (20° < 9 < 160°) 
and forward (7° < 9 < 25°) tracking detectors are used to measure charged particle trajectories, 
to reconstruct the interaction vertex and to complement the measurement of hadronic energy. 
The innermost proportional chamber CIP (9° < 9 < 171°) is used to veto charged particles for 
the identification of photons. The LAr calorimeter and inner tracking detectors are enclosed in 
a superconducting magnetic coil with a field strength of 1.16 T. From the curvature of charged 
particle trajectories in the magnetic field, the central tracking system provides transverse mo- 
mentum measurements with a resolution of ap T /Pr = 0.5% Px/GeV © 1.5% [35]. The return 
yoke of the coil is the outermost part of the central detector and is equipped with streamer tubes 
forming the central muon detector (4° < 9 < 171°). In the forward region of the detector 
(3° < 9 < 17°) a set of drift chambers detects muons and measures their momenta using an 
iron toroidal magnet. The luminosity is determined from the rate of the Bethe-Heitler process 
ep^ep^f, measured using a photon detector located close to the beam pipe at z = —103 m, in 
the backward direction. 



5 Data Analysis 

The triggers employed for collecting the data used in this analysis are based on the detection of 
electromagnetic and hadronic energy deposits or missing transverse energy in the LAr calorime- 
ter [36]. For events with missing transverse energy of 20 GeV, the trigger efficiency is ~ 90% 
and increases to above 95% for missing transverse energy above 30 GeV. Events containing an 
electromagnetic deposit (electron or photon) with an energy greater than 10 GeV are triggered 
with an efficiency close to 100% [37]. Events with two or three jets of transverse momentum 
larger than 20 GeV are triggered with an efficiency of nearly 100%. 

In order to remove background events induced by cosmic showers and other non-ep sources, 
the event vertex is required to be reconstructed within 35 cm in z of the nominal interaction 
point. In addition, topological filters and timing vetoes are applied [38]. 

In a first analysis step, calorimetric energy deposits and tracks of the event are used to look 
for electron, photon and muon candidates. Electron and photon candidates are characterised 
by compact and isolated electromagnetic showers in the LAr calorimeter. The identification of 
muon candidates is based on a track measured in the inner tracking systems associated with sig- 
nals in the muon detectors [39,40]. Calorimeter energy deposits and tracks not previously iden- 
tified as electron, photon or muon candidates are used to form combined cluster-track objects, 
from which the hadronic energy is reconstructed [41,42]. Jet candidates are reconstructed, with 
a minimum transverse momentum of 2.5 GeV, from these combined cluster-track objects using 
an inclusive k T algorithm [43,44] with a P T weighted recombination scheme in which the jets 
are treated as massless. The missing transverse momentum P™ 1SS of the event is derived from all 

2 This device was installed in 1995, replacing a lead-scintillator "sandwich" calorimeter [30]. 
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detected particles and energy deposits in the event. In events with large P™ 1SS , a neutrino can- 
didate is reconstructed. The four-vector of this neutrino candidate is calculated assuming trans- 
verse momentum conservation and the relation ^^E 1 — P*) + [E v — P") = 2E® = 55.2 GeV, 
where the sum runs over all detected particles; P z is the momentum along the proton beam axis 
and is the electron beam energy. The later relation assumes that no significant losses are 
present in the electron beam direction. 

In a second step, additional requirements are applied to ensure a clear identification of par- 
ticles. For electrons and photons, the hadronic energy within a distance in the pseudorapidity- 
azimuth (rj — 0) plane R = a/ Arj 2 + A0 2 < 0.5 around the electron (photon) is required to be 
below 3% of the electron (photon) energy. Furthermore, each electron (photon) must be isolated 
from jets by a minimum distance in pseudorapidity-azimuth of R > 0.5 to any jet axis. In the 
polar angle region 9 e > 35° electrons must be associated to a charged track and be isolated from 
any other track by a minimum distance of R > 0.5. In the central region (# 7 > 20°), photons are 
selected only if no track points to the electromagnetic cluster neither within a distance of closest 
approach (DCA) of 30 cm nor within R < 0.5. In the forward region (9 1 < 20°) only photon 
candidates with no good quality track with a DCA to the cluster below 12 cm are accepted. 
In this region, the calorimetric isolation of the photon candidate is tightened by requiring that 
the hadronic energy within R < 1, instead of R < 0.5, around the photon be below 3% of 
the photon energy. In addition, it is required that no hit in the CIP be associated to the photon 
candidate. A muon should have no more than 5 GeV deposited in a cylinder, centred on the 
muon track direction, of radius 25 cm and 50 cm in the electromagnetic and hadronic sections 
of the LAr calorimeter, respectively. Additionally, the muon is required to be separated from 
the closest jet and from any track by R > 1 and R > 0.5, respectively. 

Specific selection criteria applied in each decay channel are presented in the following sub- 
sections. 

5.1 qj Resonance Search 

The signature of the q*— decay channel consists of one high P T isolated electromagnetic 
cluster and one high P T jet. SM background arises from radiative NC DIS and prompt pho- 
ton events. As decay products of a massive particle have large transverse momenta and are 
boosted in the forward region, events are selected with a photon with transverse momentum 
Pj. > 35 GeV in a polar angle range 5° < 6 1 < 90°. The events are required to have at least 
one jet in the polar angle range 5° < 9 ]Ct < 80° with a transverse momentum larger than 
20 GeV. Photoproduction background events typically have low values of the Bjorken scaling 
variable, x h , calculated from the hadronic system using the Jacquet-Blondel method [45,46]. 
Their contribution is reduced by a factor of two by requiring x h > 0.1. Finally, to further re- 
duce the background from NC DIS, it is required that no electromagnetic deposit with an energy 
larger than 10 GeV with an associated track is present in the LAr. 

After this selection, 44 events are found in the data compared to a SM expectation of 46 ± 8 
events. The errors on the SM prediction include model and experimental systematic errors 
added in quadrature (see section 5.5). The remaining dominant SM background contributions 
are prompt photon (66%) and radiative NC DIS (26%) events. The invariant mass of the q* 
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candidate is calculated from the four-vectors of the photon and the jet candidate of highest 
P T . The invariant mass distribution of the q* candidates and the SM background expectation is 
presented in figure 2(a). From Monte Carlo studies, the selection efficiency is 35% for M q * = 
120 GeV, increasing to 45% for M q * = 260 GeV. The total width of the reconstructed q* mass 
distribution is 6 GeV for a generated q* mass of 120 GeV, increasing to 12 GeV for a q* mass 
of 260 GeV 



5.2 qqq Resonance Search 

The signatures of the q*-^qW^qqq and q*^>qZ^qqq decay channels are similar to each other 
and consist of three high transverse momentum jets. The SM background is dominated by 
multi-jet photoproduction and NC DIS events. Events are selected with at least three jets in the 
polar angle range 5° < 9 ]Ct < 120° with transverse momenta larger than 50, 30 and 15 GeV, 
respectively. In each event, a W or Z boson candidate is reconstructed from the combination of 
the two jets with an invariant mass closest to the nominal W or Z boson mass. The reconstructed 
mass of the W or Z candidate is required to be larger than 70 GeV and smaller than 100 GeV. 
From MC studies, in decays of q* of large mass, the highest Pt jet often does not originate from 
the boson decay. Therefore, only events in which the highest P T jet is not associated to the W 
or Z boson candidate are selected. This requirement is particularly effective in suppressing the 
photoproduction background at high q* masses. However, it reduces the q* selection efficiency 
at low masses. 

After this selection, 341 events are observed while 326 ± 78 are expected from the SM. The 
remaining dominant SM background contributions are photoproduction (52%) and NC DIS 
(39%) events. The invariant mass of the q* candidate is calculated from the highest P T jet and 
W or Z candidate four-vectors. The invariant mass distributions of the q* candidates and of 
the SM background are presented in figure 2(b). The selection efficiency in this channel is 5% 
for M q * = 120 GeV, increasing to 35% for M q « = 160 GeV and to 55% for M q * = 260 GeV. 
The total width of the reconstructed q* mass distribution is 11 GeV for a generated q* mass of 
120 GeV, increasing to 25 GeV for a q* mass of 260 GeV. 

5.3 qev and q^iu Resonance Searches 

The signature of the q*^qW^qeu and q*^qW^qfiu decay channels consists of one ener- 
getic lepton, a prominent jet and missing transverse momentum. Events with this topology 
correspond in the SM to single W production [47]. Other SM background processes that may 
mimic the signature through misidentification or mismeasurement are NC and CC DIS, photo- 
production and lepton pair production. 

In the search for q*^qW^qeu, events with P™ 1SS > 25 GeV, one electron with P| > 
10 GeV and one jet with P^ ct > 20 GeV are selected. The electron and the jet must be detected 
in the polar angle range 5° < 6 e < iCt < 100°. Furthermore, the electron must be isolated from 
jets by a minimum distance of R > 1. The ratio V ap /V p of transverse energy flow anti-parallel 
and parallel to the hadronic final state [45] is used to suppress photoproduction and NC DIS 
events. Events with V ap /V p > 0.25 are rejected. The invariant mass of the W boson candidate, 
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reconstructed from the four-vectors of the electron and neutrino candidates, is required to be 
between 55 and 100 GeV. After this selection six data events remain, while 6.0 ± 0.8 SM 
events are expected, of which 82% are from single W production. The invariant mass of the 
q* candidate is calculated from the jet and W candidate four-vectors. For this calculation, 
the W candidate four-vector is scaled such that its mass is set to the nominal W boson mass. 
The invariant mass distribution of the q* candidates and the SM background is presented in 
figure 2(c). The selection efficiency in this channel is ~ 20% for M q * = 120 GeV, increasing to 
30% for M q * = 260 GeV. The total width of the reconstructed q* mass distribution is 10 GeV 
for a generated q* mass of 120 GeV, increasing to 20 GeV for a q* mass of 260 GeV. 

In the search for q*^qW^qfiu, events with P™ 1SS > 25 GeV, one muon with P£ > 10 GeV 
and one jet with P^ ct > 15 GeV are selected. The muon and the jet must be detected in the polar 
angle ranges 5° < 9^ < 100° and 5° < # jct < 160°, respectively. To reduce the background 
contribution from SM W production, the P T of the jet is required to be larger than 25 GeV 
in the polar angle range 9 iCt < 60°. A W candidate is reconstructed from the combination of 
the four-vectors of the muon and neutrino candidates and its mass is required to be larger than 
40 GeV. After this selection five data events remain, while 4.4 ± 0.7 SM events are expected, of 
which 90% are from single W production. The invariant mass of the q* candidate is calculated 
from the jet and W candidate four- vectors. For this calculation, the W candidate four-vector is 
scaled such that its mass is set to the nominal W boson mass. The invariant mass distribution of 
the q* candidates and the SM background is presented in figure 2(d). The selection efficiency 
in this channel is ~ 20% for M q * = 120 GeV, increasing to 40% for M g * = 260 GeV. The total 
width of the reconstructed q* mass distribution is 14 GeV for a generated q* mass of 120 GeV, 
increasing to 30 GeV for a q* mass of 260 GeV. 

5.4 qee and qfifi Resonance Searches 

The signature of the q*^qZ—^qee and q*^qZ^qfifi decay channels consists of two high P T 
leptons and an energetic jet. The production of lepton pairs constitutes the main background 
contribution from SM processes [40]. 

In the search for q*^qZ^qee, events with two electrons and one jet of high transverse mo- 
menta are selected. Events are selected with two electrons in the polar angle range 5° < 9 e < 
100° and transverse momenta larger than 20 and 10 GeV. A jet with a transverse momentum 
larger than 20 GeV must be detected in the polar angle range 5° < 9 iCt < 100°. To reduce the 
background from QED Compton and NC DIS processes, each electron must be associated to a 
good quality track also in the forward region (5° < 9 e < 35°). A Z candidate is reconstructed 
from the combination of the two electrons and its reconstructed mass is required to be compat- 
ible with the nominal Z boson mass within 7 GeV. After this selection no data event remains, 
while 0.44 ± 0.08 SM events are expected. 

In the search for q*^qZ^q/ifi, events with two muons and one jet of high transverse mo- 
menta are selected. Events are selected with two muons in the polar angle range 5° < 9^ < 160° 
and transverse momenta larger than 15 and 10 GeV, respectively. A jet with a transverse mo- 
mentum larger than 20 GeV must be detected in the polar angle range 5° < 9 ]Ct < 100°. A Z 
candidate is reconstructed from the combination of the two muons and its reconstructed mass is 
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required to be larger than 50 GeV. After this selection no data event remains, while 0.87 ± 0.11 
SM events are expected. 

In both channels, the selection efficiency is ~ 30% for events with m q * > 120 GeV. The 
total width of the reconstructed q* mass distribution in the gee (g/i/i) channel is 5 GeV (25 GeV) 
for a generated q* mass of 120 GeV, increasing to 20 GeV (30 GeV) for a q* mass of 260 GeV. 

5.5 Systematic Uncertainties 

The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies between 0.7% and 2% de- 
pending on the polar angle [48]. The polar angle measurement uncertainty is 3 mrad 
for electromagnetic clusters. The identification efficiency of electrons is known with an 
uncertainty of 3%. 

• The efficiency to identify photons is known with a precision of 10% for photons with 
P T > 10 GeV [48]. 

• The scale uncertainty on the transverse momentum of high P T muons amounts to 2.5%. 
The uncertainty on the reconstruction of the muon polar angle is 3 mrad. The identifica- 
tion efficiency of muons is known with an uncertainty of 3%. 

• The hadronic energy scale is known within 2% [48]. The uncertainty on the jet polar 
angle determination is 10 mrad. 

• The uncertainty on the trigger efficiency is 3%. 

• The luminosity measurement has an uncertainty of 3%. 

The effect of the above systematic uncertainties on the SM expectation and the signal ef- 
ficiency is determined by varying the experimental quantities by ±1 standard deviation in the 
MC samples and propagating these variations through the whole analysis chain. 

Additional model systematic uncertainties are attributed to the SM background MC gen- 
erators described in section 3. An error of 20% is attributed to NC DIS, CC DIS and photo- 
production processes with at least two high P T jets. It includes uncertainties from the proton 
distribution functions, from missing higher order QCD corrections and from hadronisation. The 
error on the elastic and quasi-elastic QED Compton cross sections is conservatively estimated 
to be 5%; the error on the inelastic QED Compton cross section is 10%. The errors attributed 
to lepton pair and W production are 3% and 15%, respectively. The total error on the SM back- 
ground prediction is determined by adding the effects of all model and experimental systematic 
uncertainties in quadrature. 

The theoretical uncertainty on the q* production cross section is dominated by the uncer- 
tainty on the scale at which the proton parton densities are evaluated. It is estimated by varying 
this scale from V§/2 to 2\f§. The resulting uncertainty depends on the q* mass and is 5% at 
M g * = 100 GeV, increasing to 12% at M q * = 300 GeV. 
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6 Interpretation and Limits 



The event yields observed in all decay channels are in agreement with the corresponding SM 
expectations and are summarised in table 1 . The SM predictions are dominated by photoproduc- 
tion and NC DIS for searches in the q^y and qqq channels and by SM W production in the qeu 
and q/ju channels. The observed invariant mass distributions are in agreement with those of the 
SM background as shown in figure 2. No data events are observed in channels corresponding 
to leptonic decays of the Z boson, in agreement with the low SM expectations. 

Since no evidence for the production of excited quarks is observed, upper limits on the q* 
production cross section and on the model parameters are derived as a function of the mass of 
the excited quark. Limits are presented at the 95% confidence level (CL) and are obtained from 
the mass spectra using a modified frequentist approach which takes statistical and systematic 
uncertainties into account [49]. 

Upper limits on the product of the q* production cross section and of the q* final state 
branching ratio are shown in figure 3. The analysed q* decays into W and Z gauge bosons 
are combined. The resulting limit on //A after combination of all decay channels is displayed 
as a function of the q* mass in figure 4, for f s — and the conventional assumption / = /'. 
Under the assumption f /A = 1/M q * excited quarks with masses below 252 GeV are excluded. 
The individual limits from different q* decay channels are also shown in figure 4(a). At low 
mass, the combined limit on //A is dominated by the q*^q^ channel, while the q* — > qW and 
q* — > qZ channels start to contribute for masses above 150 GeV and dominate for masses above 
200 GeV. These new results extend the previously published limits by HI [2] and ZEUS [3] 
by a factor of two to five in //A. Constraints on q* masses beyond the HERA kinematic 
limit are set. Excited quarks with masses below 380 GeV are excluded for coupling values 
f/A > 0.03 GeV- 1 (see figure 4(b)). The exclusion limit obtained at LEP by the DELPHI 
Collaboration [4] assuming that the branching ratio of the q* — > q^y is equal to 1 is also shown 
in figure 4(b) and is considerably weaker than the present result. 

The sensitivity of this analysis to non-zero values of f s and its complementarity to results 
obtained at the Tevatron is also studied. The limit on / obtained for two example values of f s 
and under the assumptions f — f and A = M q * is presented in figure 5. This limit is derived 
using the 7, W and Z decay channels of excited quarks. The q* — > qg decay channel gives rise 
to a dijet resonance. It was verified in a complementary analysis [50] that no resonance from 
two high P T jets is observed in the present data. Due to the overwhelming dijet SM background, 
the total limit on q* production is not improved if the q* — > qg decay channel is included. The 
present limit is compared to the limit obtained by the CDF Collaboration for f s = 0.1 and the 
same assumptions [5]. For f s < 0.1 and for M q * < 190 GeV, the present analysis probes a 
domain not excluded by Tevatron experiments. In the case / = /' = f s = 1 and A = M q *, 
Tevatron experiments are able to exclude excited quark masses up to 870 GeV [7, 8]. 

7 Conclusion 

A search for the production of excited quarks is performed using the full e ± p data sample col- 
lected by the HI experiment at HERA with an integrated luminosity of 475 pb -1 . The excited 
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quark decay channels g*— >qj, q*^qZ and q*^qW with subsequent hadronic or leptonic de- 
cays of the W and Z bosons are considered and no indication of a signal is found. Improved 
limits on the production cross section of excited quarks are obtained. Within gauge mediated 
models, an upper limit on the coupling //A as a function of the excited quark mass is estab- 
lished for the specific relations f — f and f s = 0. For f/A = 1/M q * excited quarks with a 
mass below 252 GeV are excluded at 95% confidence level. The results presented in this paper 
extend previously excluded domains at HERA and LEP and are complementary to q* searches 
performed at the Tevatron. 
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HI Search for q* at HERA (475 pb" 1 ) 



Channel Data SM Signal Efficiency [%] 

q*^ qi 44 46 ± 4 ± 7 35-45 

q*^qW/Z^qqq 341 326 ± 48 ± 62 5-55 

g*->gW->gei/ 6 6.0 ± 0.2 ±0.8 20-30 

q*^qW^qfiv 5 4.4 ± 0.2 ±0.7 20-40 

q*^qZ^qee 0.44 ± 0.06 ± 0.04 15-30 

q*^qZ^qnn 0.87 ± 0.10 ± 0.04 15-30 



Table 1: Observed and predicted event yields for the studied q* decay channels. The first and 
second errors on the SM predictions correspond to experimental and model systematic errors, 
respectively. Typical selection efficiencies obtained from MC studies for q* masses ranging 
from 120 to 260 GeV are also indicated. 
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Search for q* at HERA (475 pb" 1 ) 
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Figure 2: Invariant mass distribution of the q* candidates in the q*^q r y (a), q*^qW/Z^qqq 
(b), q*— >qW— >qev (c), and q*^qW-^qfiu (d) search channels. The points correspond to the 
data and the histograms to the SM expectation after the final selections. The error bands on 
the SM prediction include model uncertainties and experimental systematic errors added in 
quadrature. The dashed line represents the reconstructed mass distribution of MC q* signal 
events with M q * = 240 GeV, with an arbitrary normalisation. 
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Search for q* at HERA (475 pb" 1 ) 
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Figure 3: Upper limits at 95% CL on the product of the q* cross section and decay branching 
ratio, ax BR, in the three types of final states for q* events as a function of the excited quark 
mass. The q* decay channels into the W and Z bosons are combined. Values above the curves 
are excluded. 




Figure 4: Exclusion limits at 95% CL on //A as a function of the mass of the excited quark 
with the assumptions f = f and f s — 0. The individual contributions of the q* decay channels 
are presented in (a). Values of the couplings above the curves are excluded. The excluded 
domain based on all HI e ± p data is represented in (b) by the shaded area. It is compared to the 
exclusion limit obtained at LEP by the DELPHI Collaboration [4] (dashed line), assuming that 
the branching ratio of the q* — > q^y is equal to 1. The curve //A = l/M q * is indicated in (b). 
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Figure 5: Exclusion limits at 95% CL on the coupling / as a function of the mass of the ex- 
cited quark, assuming A = M q * and f — f for two different values of f s . Values above the 
curves are excluded. Also shown is the exclusion limit obtained at the Tevatron by the CDF 
Collaboration [5] derived for f s = 0.1 (hatched area). 
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